Linear Langevin Equation of Critical Fluctuation in Chiral Phase Transition 
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We derive the linear Langevin equation that describes the behavior of critical fluctuation above the 
critical temperature of the chiral phase transition in the Nambu-Jona-Lasinio model. The Langevin 
equation relaxes exhibiting oscillation and shows thermalization. The relaxation becomes slower as 
the system approaches the critical point (the critical slowing down). The time correlation function 
also calculated using the Langevin equation shows that there exists the same soft mode as the 
previous calculation in the linear response theory. 
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PACS numbers: 05.10.Gq,11.10.Wx,11.30.Rd 

It is widely believed that hadronic matters undergoes a 
phase transition at high temperature, and hence becomes 
the quark-gluon plasma (QGP) where quarks and gluons 
are deconfined and chiral symmetry broken in the hadron 
phase is restored. In order to explore the QGP, relativis- 
tic heavy-ion collisions are a significant experimental in- 
strument. Heavy- ion collisions are essentially nonequi- 
librium processes and thus the description of the time 
evolution is inevitable to understand the phenomena in a 
comprehensive way. One of the powerful tools to describe 
time dependent processes is the hydrodynamic model. In 
some instances, the hydrodynamic model fairly correctly 
describes the experimental data of heavy- ion collisions, 
for example, transverse momentum dependence of the 
elliptic flow coefficient v% Q- 

Recently, our attention has been focused on the 
dynamics near the critical point of the chiral phase 
transition|2j. Unfortunately, the hydrodynamic model 
is sometimes inadequate to describe phenomena near the 
critical point. The hydrodynamic model is valid only 
in the case where we are interested in gross variables 
associated with macroscopic time and length scales (en- 
tropy density, energy density, etc) and they are widely 
separated from microscopic variables (degrees of freedom 
of each quark and gluon). However, near the critical 
point of the second order phase transition, the indefinite 
increase of the correlation range of critical fluctuations 
destroys the simple picture of widely separated scales. 
This situation is similar to the case of the classical dilute 
gas, where the hydrodynamic model loses its validity in 
describing physics with mesoscopic scales (for example, 
mean free path), and a kinetic equation becomes appro- 
priate. To describe the critical dynamics, it is probable 
that we still do not need the full details of the Heisenberg 
equation of motion because of the universality of phase 
transition. Then, from an analogy to gas dynamics, the 
critical dynamics will be described by a kinetic equation, 
like a Langevin equation^. The purpose of the present 
paper is to derive a linear Langevin equation to describe 
the dynamics of the critical fluctuation above the critical 
temperature of the chiral phase transition at zero density. 



As a low-energy effective model of QCD, we adopt the 
two-flavor and three-color Nambu-Jona-Lasinio (NJL) 
model 0; 

H = J d 3 xq(-iy-V)q-g J d 3 x{(qq) 2 + (qi l5 r) 2 } 



= Hq + Hj. 



(1) 



The coupling constant g and the three momentum cut- 
off A are fixed so as to reproduce the pion decay con- 
stant /,r = 93 MeV and the chiral condensate (qq) = 
(-250 MeV) 3 : g = 5.01 GeV~ 2 and A = 650 MeV. The 
critical temperature of the chiral phase transition using 
this parameter set is T c = 185 MeV in the mean field 
approximation. 

In principle, the dynamics of the chiral phase tran- 
sition is described by solving the Heisenberg equation 
of motion for the order parameter, that is, the chiral 
condensate (a) = (qq). However, the Heisenberg equa- 
tion incorporates motion on microscopic and macroscopic 
scales on an equal footing. To carry out coarse-graining 
of microscopic variables systematically, we introduce pro- 
jection operators P and Q = 1 - P [1 H 0, & 0. This 
technique was originally introduced by NakajimaQ to 
derive a master equation. The application to a Langevin 
equation was attempted by Mori[6j. Recently, the uni- 
fication and generalization of both treatments has been 
realized 0> @- It should be noted that the elimination of 
microscopic variables does not mean neglecting their ef- 
fect. They give rise to dissipation terms and noise terms 
in the coarse-grained macroscopic equations. 

Then, the time evolution of the chiral phase transition 
is described by @ 



dt 



Sa{x,t) = e lLt PiL5a(x,0) 



f dse lL{t - s) PiLQe lQLoQs iL6a(x,0)+£(^tp) 
Jo 



where L and Lq are Liouville operators of the Hamil- 
tonian H and Hq, respectively. Here, we calculate the 
fluctuation of the order parameter from its equilibrium 
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value: <5><t(x, t) = cr(x, t)— < cr(x) > eq . This equation is 
called the time-convolution (TC) equation0,Q- The first 
term on the r.h.s of the equation expresses the term corre- 
sponding to a collective oscillation such as plasma wave, 
spin wave, etc. The second term is the memory term 
that causes dissipation. Here, we have already expanded 
the memory term up to lowest order in Hj, following 
Ref. ||. The third term is the noise term and we do not 
give its concrete form here. This is defined through the 
fluctuation-dissipation theorem of second kind (2nd F-D 
theorem) as we will see later. 

There are several possible projection operators 
that extract slowly varying parts from an operator 
H H E3, E3 In this paper, we adopt the Mori 
projection operator(MPO) that projects any operators 
onto the space spanned by gross variables j6j. To define 
it, we must realize the time scale of gross variables 
are extremely long compared to that of other micro- 
scopic variables. There are three candidates for gross 
variables^: (i) order parameters, (ii) density variables 
associated with conserved quantities and (iii) their 
products. However, near critical points of the second 
order phase transition, the order parameters are proba- 
bly much slower than the density variables because of 
critical slowing down(CSD). In this sense, we can ignore 
the time evolution of density variables and exclude (ii) 
from our gross variables. Furthermore, we ignore (iii) 
for simplicity. We shall discuss the importance of (iii) 
at the end of this paper. After all, the gross variable 
relevant in our calculation is the order parameter of the 
chiral phase transition tr(x, t). Then, the MPO in our 
calculation is defined by 



PO 



= J d 3 xd 3 x'{0,5a{x)) • (^(x),^'))- 1 •<5a(x'),(3) 

where O is an arbitrary operator. The inner product 
means the canonical correlation, 



X ,Y)= ^ ^-Tr{pe XH °Xe- XH °Y], 
Jo P 



(4) 



where p = exp(-/3 J ff )/Tr[exp(-/3 J ff )] and = l/T. 
The system attains the thermal equilibrium state with 
this temperature T. In the classical limit, the canonical 
correlation coincides with the classical correlation func- 
tion Tr [pXY]. 

Substituting the MPO into the TC equation ©, we 
have 



d 



6a(k,t) = - [ drr(k,r)<5a(k,t-r) +£(M). (5) 
Jo 



Here, the memory function r(k, t) is given by the inverse 
Laplace transformation of T L (k, s) , defined by 

r L (k, s) = - S ^(k)/( Xo (k) + *f(k)) • (1 - 2 ff /3 X o(k)),(6) 




FIG. 1: The time evolution of T(0,t) in the unit of A 2 . 
The solid, dashed and dotted lines are the r(0,t)/A 2 for the 
temperatures e — (T—T c )/T c = 0.05, 0.1 and 0.2, respectively. 



where the index L means the Laplace transformation, 
Xj(k) = /" dte- st d X t(k)/dt with 



Xt(k) = 



V {(1 - 1n(E )) ^p-Ep+k + P(P + k ) 



i(E p +E p+k )t , -i(_E p +£ p+k )n 



x (e 



+M£p +k ) - n(Ep)) 



ffpEp+k - p(p + k) 
EpE p+ k(Ep — -Ep+k) 



Here, n(E) = [exp((3E) + l]" 1 , E w = |k|, A^ c = 3 and 
Nf — 2. It should be noted that there is no chiral con- 
densate above T c . In the calculation of the memory term, 
we calculated only the term corresponding to a ring di- 
agram in the scalar channel and employed the technique 
proposed in Ref. [l2| . 

The memory function r L (0, s) vanishes at T c . This is 
shown by using the self-consistency condition of the chiral 
condensate calculated in the mean field approximation, 



l-2fl8xo(0)| 



T=T C 



0. 



where T c = 185 MeV. From Eqs. © and ©, one can see 
that r L (0, s) vanishes at T c . The small r L (0, s) near T c 
means the slowing down of the relaxation, that is nothing 
but CSD. If we take account of the terms beyond the ring 
diagram approximation in the calculation of the memory 
function, the critical temperature deviates from that in 
the mean field approximation. 

The TC equation (JSJ is still inadequate to be inter- 
preted as a Langevin equation because the memory func- 
tion r(k, t) converges to a finite value at t — > oo, as 
is shown in Fig. 1. As a matter of fact, it is usually 
expected that the memory function tends to zero with 
the shorter time scale than the macroscopic one, because 
we can prove that the memory function is given by the 
time correlation of the noise. This is called the 2nd F- 
D theoremQ. The anomalous behavior of the memory 
function gives rise to the oscillation in the time evolution 
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of 5cr(k,i), as we will see later in Fig. 2. This anomaly 
may arise because the MPO defined in Eq. J3J is incom- 
plete such that the noise still incorporates gross variables 
other than Sa(x). As an example, let us consider the 
Langevin equation with two gross variables x and p.|l3l| 

j f x(t) = pit), (9) 
d /"* 

-p(t) = -io 2 x{t) - J dsE(t - s)p(s) + f(t), (10) 

where is the memory function given by the time 
correlation of the noise f(t) and has no long time corre- 
lation. Solving the first equation and substituting into 
the second one, we have 

j t p(t) = - £ ds{u 2 + E(t - s))p(8) + f{t), (11) 

where we set the initial condition x(0) = 0. Although 
the equation has only one gross variable p, the memory 
function uj 2 + E(t) has the long time correlation due to 
the oscillation. In this case, the true memory term given 
by the 2nd F-D theorem is not uj 2 + E(t) but E(t). 

Thus, to define a renormalized memory function with- 
out long time correlation, we separate the oscillation ef- 
fect from the memory function. By carrying out the 
Fourier transformation, one can easily recognize that 
the correction to the frequency shift and the damp- 
ing are given by the imaginary part and the real 
part of the Fourier transform of the memory function, 
respectively Qj. Thus, the Langevin equation with the 
renormalized memory function <I>(k, i) is expressed as, 

rl 

— Sa(k,t) = dTnl{t~T)5a(k,T) 
dt J 




t 

dr$(k,r)5o-(k,t-r) + ^(k,t), (12) 



where Q k (t) = i J dujlm[T L (k, -iu> + e)]e~ iart /2tt and 
$(M) = JdcjRe[r L (k, -iuj + e)]e- iut /2tt. This mem- 
ory function does not have a long time correlation any 
more. Here, we artificially removed the oscillation effect 
from the memory term. However, a similar procedure 
should be automatically implemented by using the MPO 
defined by a complete set of gross variables. 

Now, we can define the noise with short time correla- 
tion. The definition of the noise and the 2nd F-D theorem 
lead to the following correlation properties: 

(£(k,i),fc(k',0)) = (£(k,i)> = 0, (13) 

(£(MUt (k ',f')) = vsl 3 lMk,t-t') Xo (k), (14) 

where (O) = Tr[p O] and V is the volume of the system. 
Here, we assume that the noise has the translational in- 
variance in space and time. The first correlation indi- 
cates that the noise does not include components that 
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FIG. 2: The time evolution of Sa(0, t). The solid, dashed and 
dotted lines are the Sa(0, t) for the temperatures e = 0.05, 0.2 
and 0.5, respectively. 



vary with the gross time scale and it is possible to regard 
it as a random field for 5<r(k, t). The second correla- 
tion characterizes the noise as a random field. In the 
following, we solve Eq. I|12fl as a classical equation with a 
random noise. For this, we introduce the classical noise 
that reproduces the correlation properties defined above; 

<£(M)&r(k',0) > = <£(k,t)»=0, (15) 
<£(M)r(k',t')> = (£(MU f (k',0), (16) 

where <C 3> means the average for noise with a suit- 
able stochastic weight. It is worth emphasising that the 
correlations determined here coincide with the condition 
of themalization for a Langevin equation^). In this 
sense, the system described by the Langevin equation 
approaches a thermal equilibrium state with time. 

The averaged time evolution of the critical fluctuation 
at vanishing momentum is shown in Fig. 2. One can 
see that the nonequilibrium fluctuation relaxes with os- 
cillation and finally converges zero. This indicates that 
the critical dynamics of the chiral transition may not 
be described by a diffusion equation like the time de- 
pendent Ginzburg-Landau (TDGL) equation. The rate 
of the relaxation becomes slower as the temperature ap- 
proaches T c because of CSD. The relaxation time is char- 
acterized by r rt = 2/70 where 27k = J Q ds<I>(k, s). At 
e = (T - T c )/T c = 0.2, r rt is about 10 fm, that is the 
same order as the expected life time of QGP. Thus, we 
cannot ignore the fluctuation of the order parameter at 
a temperature lower than 222 MeV. On the other hand, 
at higher temperature, the fluctuation relaxes firstly and 
other gross variables become important. 

As we have shown so far, the Langevin equation reveals 
CSD and thcrmalization. This means that the Langevin 
equation correctly describes the behavior near T c . At 
the same time, as was pointed out in Ref. 0], a soft 
mode appears above T c in the correlation function of the 
fluctuation of the order parameter. To investigate the 
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FIG. 3: The temperature dependence of 1(0, w) in the unit 
of A 2 . The solid, dashed and dotted lines are the I(0,u)/A 2 
for the temperatures e = 0.05, 0.2 and 0.5, respectively. 



there exists a soft mode. As a result, we can conclude 
that the Langevin equation fulfils the requirements near 
T c : CSD, thermalization and soft mode. 

We have discussed here only the linear Langevin equa- 
tion and ignored the nonlinear effect. However, the mode 
coupling theory makes it clear that near critical points 
the nonlinear effect becomes important because of the 
large correlation length and leads to the deviation from 
the van Hove theory in calculating dynamical critical 
exponent 0. In order to take nonlinear terms into ac- 
count, we must choose them as gross variables in defining 
the MPO. Another intriguing subject is to apply this for- 
mulation to finite density and the color superconducting 
phase transition. These are future projects. 



correlation function, we define the power spectrum as 

I(k,u>)= lim -^«|fc(k,^)| 2 », (17) 
where T is the time when we observe the system, and 

IWlr , ,M 2 ^ TV ^ [^ 2 $(k,^)xo(k)] 

< \da{k,u})\ > = IV- 2 - rr£18) 

| — uj z + 12£ — iu)<P(k, uj)\ z 

The Wiener-Khinchin theorem tells us that the correla- 
tion function is given by the power spectrum 

C(x,i) = lim < Sa(x + x', t + t')Sa(x', t') > 

t' — ^OO 



7^ /(k '" )e 



— iujt gik,x 



(19) 



The temperature dependence of the power spectrum at 
vanishing momentum is shown in Fig. 3. We can see 
that the peak moving toward origin becomes prominent 
as the temperature is lowered toward T c . The power spec- 
trum characterizes the space-time correlation in energy- 
momentum space and hence can be interpreted as the 
spectral function in the thermal Green's function. Then, 
the peak with narrow width reveals the existence of a col- 
lective mode whose energy tends to vanish as the temper- 
ature approaches T c . Such a mode is called a soft mode. 
The temperature dependence of the power spectrum is 
consistent with the previous result, where the spectral 
function is calculated in the linear response theory 15]. 

We have derived the linear Langevin equation that de- 
scribes the dynamics of the chiral phase transition in the 
NJL model using the projection operator method. The 
equation reveals CSD and shows thermalization. The or- 
der parameter relaxes exhibiting oscillation. This means 
that a simple diffusion-type equation like the TDGL 
equation may be inadequate to describe the dynamics of 
the chiral transition. The power spectrum was also cal- 
culated using the Langevin equation and we found that 
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